Temperature dependence of the nitrogen-vacancy magnetic resonance in diamond 
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The temperature dependence of the magnetic resonance spectra of nitrogen- vacancy (NV") 
ensembles in the range of 280-330 K was studied. Four samples prepared under different 
conditions were analyzed with NV" concentrations ranging from 10 ppb to 15 ppm. For all 
samples, the axial zero-field splitting (ZFS) parameter, D, was found to vary significantly 
with temperature, T, as dD/dT = —74.2(7) kHz/K. The transverse ZFS parameter, E, was 
non-zero (between 4 and 11 MHz) in all samples, and exhibited a temperature dependence 
of dE/(EdT) — —1.4(3) x 10~ 4 K _1 . The results might be accounted for by considering 
local thermal expansion. The temperature dependence of the ZFS parameters presents a 
significant challenge for diamond magnetometers and may ultimately limit their bandwidth 
and sensitivity. 



Magnetometers based on nitrogen-vacancy (NV) 
ensembles in diamond [IrEI promise high-sensitivity, 
rivaling those of superconducting quantum interfer- 
ence devices (SQUIDs) |U and alkali vapor magne- 
tometers |5], in a scaleable solid state system that 
can be operated over a wide range of temperatures. 
This remarkable combination of spatial resolution 
[HJ [7] and magnetic sensitivity [5] make diamond 
magnetometers promising candidates for remote- 
detection and low-field nuclear magnetic resonance 
spectroscopy [9Hl2j. nano-scale biological imaging 
[6l [1"3TH5| , and studies of novel magnetic and super- 
conducting materials [3j [16]. Until now, the temper- 
ature dependence of the magnetic resonance spectra 
has not been systematically studied and has only 
briefly been mentioned in the literature [3l [T7]- In 
this Letter, we report a striking temperature depen- 
dence of the magnetic-resonance spectra of NV" en- 
sembles in diamond over the temperature range of 
280-330 K. These findings have important implica- 
tions for the design of diamond magnetometers and 
may ultimately limit their sensitivity. 

The resonance spectra were recorded using the 
continuous-wave Fluorescence Detected Magnetic 
Resonance (FDMR) method Q71 [IB]. Light from 
a 514-nm Argon-ion laser was focused with a 2.5 cm 
focal length lens onto the diamond samples, excit- 
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ing the NV" centers' 3 ^2 — > optical transition 
via a phonon sideband |19j . The same lens was 
used to collect fluorescence from the diamond which 
was then passed through a dichroic mirror and a 
650-800 nm bandpass filter and detected with a pho- 
todiode. Noise due to laser power fluctuations was 
reduced by normalizing the fluorescence signal to 
a reference photodiode which monitored the inci- 
dent laser power. The output of a microwave signal 
generator was amplified, passed through a straight 
~ 200- /im diameter copper wire of length ~ 5 mm 
placed within 500 /im of the focused light beam, 
and terminated with 50-f2 impedance. For temper- 
ature control, the diamond was thermally connected 
to a copper heatsink and placed inside an insulated 
aluminum housing. The temperatures of the heat 
sink and housing were controlled with separate ther- 
moelectric (TE) elements. Unless otherwise stated, 
the results reported in this Letter were obtained 
with a magnetic field of < 1 G, laser-light power of 
~ 150 mW, and microwave power (after the wire) 
of ~ 10 dBm. 

For temperature scans, the temperature of the 
copper plate in direct thermal contact with the di- 
amond was monitored with an AD590 sensor. The 
FDMR spectra were recorded with the temperature 
stabilized so that temperature excursions were less 
than 0.05 K over 5 min. In order to avoid stray mag- 
netic fields when recording the spectrum, the cur- 
rents supplied to both TE elements were chopped 
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at a frequency of 2 Hz using photoMOS circuits, 
and the spectra were recorded only when the TE 
currents were off. The process was repeated un- 
til the temperature had been scanned through the 
280-330 K range several times in both directions. 

The NV-ensemble magnetic-resonance spec- 
troscopy has been described, for example, in Refs. 
[21 [181 [20H22| ) and is only briefly summarized here. 
Optical pumping via a spin-selective decay path col- 
lects NV centers (total spin S = 1) in the \m s = 0) 
ground-state magnetic sublevel [19]. In the absence 
of external fields, the \m s = 0) and \m s = ±1) levels 
are split by an energy equal to the axial zero-field 
splitting (ZFS) parameter, D ~ 2.87 GHz. For per- 
fect C^ v symmetry, the transverse ZFS parameter 
is E = and the \m s = ±1) levels remain degener- 
ate. When the frequency of a microwave field that is 
transverse to the symmetry axis is tuned to the en- 
ergy splitting between the \m s = 0) and \m s = ±1) 
levels, NV centers are transferred to the \m s = ±1) 
sublevels, resulting in diminished fluorescence with 
a contrast as high as 30% p]. In the presence of 
an applied magnetic field, B, the \m s = ±1) levels 
split, revealing resonances separated by 2g^y[x^B , 
where g NV = 2.003 is the NV" Lande factor [201 125] 
and [ib is the Bohr magneton. For ensembles, there 
are four different NV orientations and, provided that 
9nv^b\B\ <C D, only the projection of the magnetic 
field on the N-V axis affects the transition frequen- 
cies [22J. 

The zero-field Hamiltonian for the ground state, 
including hyperfine coupling to the 14 N nucleus 
(spin 1=1), can be written as: 
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where A\\ = -2.1 MHz and A± = -2.7 MHz are, 
respectively, the axial and transverse hyperfine con- 
stants |23| . Analysis of this Hamiltonian reveals six 
allowed microwave transitions for each N-V orienta- 
tion. The relative intensities can be calculated by 
treating the interaction with the microwave field, 
B\ as a perturbation, M\ = qnv^bBi ■ S, with ma- 
trix elements that depend on the alignment of the 
microwave radiation with respect to the symmetry 
and strain axes of each N-V center. However, since 
the exact geometry and the number of NV" centers 
of each orientation were not known a priori, Gaus- 
sian functions with variable amplitudes and equal 
widths, centered about these transition frequencies, 
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FIG. 1: Zero-field FDMR spectrum at 293 K for S3 
and the corresponding fit based on Eq. [l] (solid green 
line). The six blue lines represent the fitted amplitudes 
at each transition frequency, and the fitted linewidth was 
3.3 MHz (full width at half maximum). The microwave 
power was reduced to ~ — 10 dBm to resolve the hyper- 
fine structure, resulting in the relatively small contrast 
of ~ 0.6%. The best-fit parameters for this scan are 
E = 4.1(2) MHz and D = 2866.8(2) MHz. 



were fit to the spectra. Including residual magnetic 
fields, measured by a commercial fluxgate magne- 
tometer to be less than 1 G, into the model did not 
significantly influence the fits. 

Four single-crystal samples of mm-scale dimen- 
sions were studied, which were labeled S2, S3, S5, 
and S8 and characterized in Ref. [2]. Figure [T] shows 
the FDMR spectrum at 293 K for S3, a sample syn- 
thesized by chemical vapor deposition (CVD) with 
[NV"]~ ~ 10 ppb [2]. As there was no applied mag- 
netic field, the splitting between resonance peaks is 
due to non-zero E, induced by local strain [6, 22, 24j. 
This feature is present in varying magnitudes for 
all four samples. Even though all four NV orienta- 
tions are present, the spectra are reasonably well- 
described by just six broad transitions, suggesting 
that the strain splittings are spatially inhomogenous 
|17[ 125] . As no correlation with NV" concentration 
was observed (see Tab. |T|, further work is necessary 
to determine the exact strain mechanism. 

During each temperature scan, the spectrum was 
fit to an empirical function similar to the one de- 
scribed above, and the ZFS parameters were ex- 
tracted. Figure [2|a) displays the spectra at two dif- 
ferent temperatures for another sample, S8, a high- 
pressure, high-temperature (HPHT) synthesized di- 
amond with [NV"] ~ 0.3 ppm, as well as the empir- 
ical fits based on Eq. [I] Fi gures ^tb) and (c) show 
the ZFS parameters as a function of temperature for 




2.85 2.86 2.87 2.88 2.89 
Microwave frequency [GHz] 









b) 



# [NV-](ppm) 



^ (kHz/K) g(MHg) ^^( IQ-^K- 1 ) 



280 290 300 310 320 330 
Temperature [K] 



2867 
2866 
2865 
2864 




280 290 300 310 320 330 
Temperature [K] 



20 


^-20 
1-40 
-60 
-80 









* dl /df f 




dD/(JT . 




* ■ 







10 20 30 40 50 60 
Laser intensity [kW/cm 2 ] 



FIG. 2: (a) Zero-field FDMR spectra at 283 K and 326 K 
for S8 with fits (solid red lines), (b) Value of E for S8 
as a function of temperature with linear fit (solid black 
line), (c) D for S8 vs. temperature with linear fit. (d) 
dD/dT and dE/dT function of laser intensity for 
S5. The dotted lines are the laser-intensity-independent 
values used in Tab. HI 

this sample. Linear least-squares fits yield dE/dT = 
-0.4(2) kHz/K and dD/dT = -76(1) kHz/K. Fig- 
ure f2Tci) displays the laser-intensity dependence for 
S5, an HPHT diamond with [NV~] ~ 12 ppm. Lin- 
ear fits (not shown) determined that any depen- 
dence of dD/dT or dE/dT on laser intensity is not 
statistically significant. Additional tests for depen- 
dence on microwave power, external magnetic field, 
and sample positioning also did not show statisti- 
cally significant effects. 

A similar procedure was performed for the three 
other samples: S2, an HPHT diamond with [NV~] ~ 
16 ppm, as well as S3 and S5 (already mentioned). 
Table |T] displays the temperature dependence of the 
ZFS parameters for each of these samples. The tem- 
perature dependence of D is similar for each sample, 
indicating that the mechanism responsible for this 
temperature variation is intrinsic to the NV cen- 
ters themselves. Taking a weighted average over 
all samples gives dD/dT = —74.2(7) kHz/K and, 
using the fitted room-temperature values for each 
sample {D ~ 2867(1) MHz), this corresponds to a 
fractional temperature dependence of dD / (DdT) = 
— 2.59(2) x 10~ 5 K^ 1 . The weighted average over 
samples of the fractional variation of E with tem- 
perature (final column of Tab. |T| is also statistically 
significant, dE/{EdT) = -1.4(3) x 1(T 4 K _1 , but 
further work is necessary to understand the nature 
of E. 

The origin of D is expected to be predominately 
due to dipolar spin-spin coupling between the two 
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TABLE I: ZFS parameters and uncertainties for four 
different samples. The values of E represent the ex- 
pected value of E(293 K) extrapolated from the lin- 
ear fits, and the error bars represent the standard error 
from the fit but not systematic effects due to imperfect 
assumptions in the model (see text). The laser inten- 
sity was ~ 25-50 kW/cm 2 throughout the collection vol- 
ume. Note that for the S2 spectra a magnetic field of 
B± ss 13 G was applied. This field enabled the isolation 
of a single NV orientation, and the simplified spectrum 
was used to verify the robustness of the model. 



unpaired electrons forming the center |18[ 120) 126] . 
This suggests a a likely mechanism for the temper- 
ature variation is local lattice expansion. Assum- 
ing that the angular electronic wavefunctions are 
temperature-independent and that D is entirely due 
to dipolar coupling, the effect of lattice expansion on 
D is: 



1 dD 
D~dT 



1 d((r 2 12 -3z 2 12 )/rl 2 )dR 



D 



dR 



dT : 



(2) 



where ryi is the displacement between the two spins, 
Zyi is the component of ryi along the N-V symme- 
try axis, and R is the distance between two basal 
carbon nuclei. The effect of thermal expansion on 
r\ 2 ) can be estimated by treating spins, 
localized near the basal carbon atoms |18| l2"T] , with 
p-orbitals |18| [23] oriented along axes 110° apart 
[28], and calculating the integral for neighboring 
values of R. Using the room-temperature values 
for bulk diamond of R = 0.252 nm and dR/dT = 
2.52 x 10~ 5 nm/K [29J, we calculate D = 2.66 GHz, 
which is within 10% of the experimental value, and 
dD/(DdT) = -5.8 x 10~ 6 K -1 , which is about a fac- 
tor of 4.5 smaller than the experimental value from 
this work. The latter discrepancy suggests that the 
macroscopic thermal expansion is not a good de- 
scription of dR/dT in the immediate vicinity of the 
defect. Ab initio calculations [26H2"8] [30] [31] which 
include the determination of local thermal expan- 
sion effects would give a more accurate prediction 
of dD/dT. 

The sharp temperature dependence of D presents 
a technical challenge for room-temperature diamond 
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magnetometry. Even if the ambient temperature 
can be controlled at the 1-mK level, this would lead 
to fluctuations in the resonance frequency of 80 Hz 
corresponding to a magnetic-field variation of 3 nT. 
Monitoring both of the Am s = ±1 resonances could 
provide a feedback mechanism for controlling this 
effect for slow drifts, since the energy difference be- 
tween these resonances does not depend on D. 

Higher-frequency temperature fluctuations due 
to, for example, laser-intensity noise, present an ad- 
ditional complication for magnetometry in the high- 
density limit. Consider the case of a Ramsey- type 
magnetometer making use of repeated light pulses 
[H El El [15] which transfer an energy to the di- 
amond on the order of E p ~ Ae[NV~]F, where 
Ae ~ 0.6 eV is the difference in energy between 
absorbed and radiated photons, V is the effective 
volume being heated, and we have conservatively 
neglected non-radiative transfer from the NV~ sin- 
glet decay path |32J and other impurities |33| . If the 
pulses are separated in time by a precession win- 
dow, r, then in steady state the diamond temper- 

, , , dT E p 

ature is modulated at a rate — — , where 

dt Vct 

c = 1.8 J/cm 3 /K is the volumetric specific heat 
of diamond |34] , Integration over the precession 
window yields a magnetometer offset of B Q ff ~ 
vrAe[NV-] dD nn m 

— w —80 nl at room temperature for 

9NVVBC dT 

[NV~] = 1 ppm. This offset makes the magne- 



tometer sensitive to laser-pulse fluctuations. Uncor- 
related, normally-distributed fluctuations in E p by 
a fraction \ produce magnetic field noise-per-unit- 
bandwidth at the level of x\B ff \j\fr ~ 1 pT/y/tlz, 
using x = 0.01 an( i t = 1 fjts. We note that this 
magnetometer noise is directly correlated with laser- 
intensity noise and therefore monitoring the incident 
laser intensity could significantly reduce this effect. 

In this work, we have measured the tempera- 
ture dependence of the ZFS parameters of four di- 
amond samples covering a wide range of NV~ con- 
centrations. We have found a significant variation 
of the axial ZFS, D, with temperature and surmise 
that it is due to local thermal expansion. We also 
present evidence of a non-zero transverse ZFS, E, 
and measure a small fractional temperature depen- 
dence just above the experimental uncertainty. The 
results have a major impact on the performance of 
NV-ensemble magnetometers and may ultimately 
limit their sensitivity and bandwidth. We expect 
that proper feedback mechanisms, such as monitor- 
ing laser intensity fluctuations and observing both 
Am s = 1 coherences simultaneously, will help to 
partially mitigate these effects. 
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